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Abstract.  The  derivation  of  moment  equations  for  the  description  of  a  flowing  binary  non-reacting  gas 
mixture  is  presented.  The  equations  obtained  are  based  on  a  kinetic  model  in  a  relaxation  form  and  can 
be  regarded  as  the  extension  of  the  quasigasdynamic  equations  (QGD)  studied  before.  As  an  example  of 
implementation  the  results  for  the  shock  wave  structure  in  an  argon-xenon  mixture  and  for  the  binary 
diffusion  process  of  argon  and  helium  are  presented  in  comparison  with  the  DSMC  results. 


I  INTRODUCTION 

There  are  two  groups  of  models  for  calculating  gas-mixture  flows.  The  first  one  consists  of  kinetic  models,  i.e. 
models  based  on  direct  numerical  simulation  methods  (DSMC)  or  on  the  solution  of  the  Boltzmann  equation 
(e.g.  [1]).  The  other  group  consists  of  systems  of  macroscopic  equations  which  are  derived  on  the  basis  of 
Navier-Stokes  equations,  in  general  by  a  phenomenological  way  (e.g.  [2]  -  [4]). 

In  this  paper  we  propose  a  new  macroscopic  model  to  describe  flows  of  a  binary  non-reacting  gas  mixture.  The 
model  is  a  two-fluid  approximation  that  consists  of  a  system  of  equations  for  density,  momentum  and  energy 
of  each  component.  The  system  of  macroscopic  equations  (named  QGDM)  is  based  on  the  kinetic  equation  in 
its  relaxation  approximation  and  is  a  natural  generalization  of  the  recently  proposed  quasigasdynamic  (QGD) 
equations  (e.g.  [5]  -  [7])  for  a  gas  mixture.  Note  that  previously  the  QGD  equations  were  generalized  for  gas 
flows  in  translational  [8]  and  rotational  [9]  non-equilibrium. 


II  KINETIC  MODEL  AND  MACROSCOPIC  EQUATIONS 

We  consider  a  mixture  consisting  of  gas  a  and  gas  b.  Each  gas  is  characterized  by  its  density  pi ,  temperature 
Tj  and  macroscopic  velocity  u,,  where  i  =  a,b.  We  introduce  the  gas  constant  i?,-  =  k/rrii,  where  k  is  the 
Boltzmann  constant,  rrii  is  the  mass  of  the  molecule.  According  to  [10]  the  kinetical  model  for  the  mixture  can 
be  expressed  as 


^  +  (^V)/a  =  ^aa(Fa  -  fa)  +  l^ab{Fa  -  fa),  (1) 

M  +  (eV)/5  =  VbbiFb  -  fb)  +  l^baiFb  -  fb),  (2) 

where  fa  and  fb  are  the  distribution  functions  for  species,  ^  is  the  molecule  velocity,  Vaa  and  Vbb  are  the 
frequencies  of  self-collisions,  Vab  is  the  frequency  of  collisions  of  a  molecules  with  b  molecules,  Vba  is  the  frequency 
of  collisions  of  b  molecules  with  a  molecules.  Fa,  Fb  and  Fa,  Fb  are  Maxwellian  distribution  functions  based 
on  [ua,Ta),  {ub,Tb),  {ua,Ta),  {ub,Tb).  The  so-called  free  parameters  (overlined  symbols)  can  be  expressed  in 
terms  of  parameters  m,u,T  of  each  gas 
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We  suppose  that  the  distribution  functions  can  be  approximated  by  a  gradient  expansions  as  in  [8] 

fi  =  Fi-  T{^V)Fi 


Ta=Ta  + 

n  =  Tb  + 


2mamb  (  mb  ^ 

-  Ta  +  —  (Ub  -  Ua) 

I  ma  ,  ^ 

[Ta  “  “  “«) 


{ma  +  mbY 
‘2mamb 

{ma  +  mbY 


(3) 

(4) 


where  r  is  identified  with  )i/p,  is  the  viscosity  and  p  the  pressure  of  the  mixture.  These  expressions  are 
substituted  for  fa  and  fb  into  the  convective  terms  of  Eqs.1-2  to  yield  an  approximation  of  these  equations 


^  +  i^V)Fa  -  {fV)T{YV)Fa 


^aa{Fa  fa}  d“  ^ab{Fa  fa}j 


(5) 


+  {f^)Fb  -  (CV)r(fV)Ff,  =  Ubb{Fb  -  fb)  +  naiFb  -  fb)-  (6) 

Equations  5  and  6  are  integrated  over  the  velocity  space  after  multiplications  by  1,  and  f^/2,  successively 
and  result  in  a  system  (QGDM)  of  macroscopic  equations  for  the  binary  mixture 

^Pa  +  ViPa  <  -  (VjPa  «  +  V‘pa)  =  0,  (7) 


d 

dt 


Pa  <  +  ViPa  uiu^a  +  =  Vi  T  {VjPa  V^aK^  +  ^  Pa  4  +  V"'  Pa  <)  +  V''  T  Vi  Pa  K  +  5“  , 


(8) 


d  ■  •  •  1  • 

-^Ea  +  ViU\{Ea  +Pa)  =  Vi  T  {V j{Ea  +  2pa)u'‘aUi  +  -V'UakUaPa)  + 


J±^Vi  T  Pa  +  Prf^^^ViTPa  V*  ^  +  5f , 

-1  Pa  7a  -  1  Pa 


iPa 


la 


Pa 


Vi  r  — ' 

7a  -  1  Pa 

where  the  volumic  energy  writes 

Ea  =  {paul)/2+Pa/{'ya  -  1)  with  Pa  =  PaRaTa- 


(9) 

(10) 


Similar  equations  hold  for  gas  b.  The  Prandtl  number  Pr  and  specific  heat  ratio  7  do  not  appear  in  the  above- 
described  treatment  and  have  been  introduced  artificially  for  generalization  purpose  to  make  these  equations 
coincide  with  Navier-Stokes  (NS)  equations  in  the  case  of  a  one-species  gas  and  vanishing  Knudsen  numbers 
(r^0)[7]. 

The  exchange  terms  S  result  from  the  integration  of  the  right-hand  sides  of  the  equations  5,  6  and  are 
recognized  as 

SI  =  VabPa{Ua  -  Ua),  Sb  =  VbaPbiub  -  Ub) ,  Sa  =  l^ab{Ea  -  Ea) ,  Sb  =  Vba{Eb  -  Eb) ,  (11) 

where  Ea  and  Eb  have  expressions  similar  to  Ea  in  Eq.lO. 

5:  +  5,“  =  0,  5f  +  5f=0  (12) 


which  is  consistent  with  momentum  and  energy  conservation. 

Thus  the  binary  mixture  has  been  described  by  a  self-contained  set  of  two-fluid  equations  which  allows  for 
different  velocities  and  temperatures  of  the  components. 

The  mixture  parameters  (no  subscript)  relate  to  the  species  parameters: 


n  =  na  +  nb,  p  =  Pa+Pb,  P  =  Pa+Pb,  U  =  {paUa  +  PbUb)/p, 

T  =  {uaTa  nbTb)/n,  m  =  {maUa  -b  mbnb)/n,  p  =  pRT, 

R=  {paRa  +  PbRb)/p  =  k/m.  (13) 
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Ill  DETERMINATION  OF  COLLISION  FREQUENCIES 


To  close  the  system  of  QGDM  equations,  it  is  necessary  to  estimate  the  cross-collision  frequencies  i^ab  and 
i^ba  and  the  relaxation  time  r. 

Relations  between  frequencies  of  collisions  of  a  molecules  with  each  other  (self-collisions)  and  with  h  molecules 
(cross-collisions),  according  to  [10],  can  be  calculated  in  the  following  way: 


^ab  — 


nig  +  mb  Hg 
2mgmb  Ub  ’ 


(14) 


where  dg  is  the  effective  molecular  diameter  for  gas  a,  dgt  is  the  effective  diameter  which  can  be  determined, 
for  example,  according  to  [1],  p.l6,  as  dgb  =  0.5(£/a  -|-  db).  In  turn,  the  collision  frequency  i^g  can  be  related  to 
with  the  gas  viscosity  in  the  approximation  of  the  VHS  and  VSS  models  for  particle  interactions  (  [1],  p.90  ): 


Pa  I  \ 

Vg  =  - n(UJg,ag), 

Pa 

where  51(wo,ao)  = 


Pa  —  Paref  (rp  )  “  ) 
aref 

5(Q;a  +  l)(ag  d-  2) 

ag{7  -  2lj<,)(5  -  2uig) ' 


(15) 


For  further  calculations  ««  =  1  is  used,  that  conforms  to  the  VHS  model  (  [1],  p.41).  In  this  case  we  write 

1)  =  ft(Wa). 

The  total  number  of  collisions  between  molecules  of  gases  a  and  b  should  be  balanced,  i.e.  relation  ng  Vgb  = 
Ub  Vba  must  be  satisfied. 

To  estimate  the  relaxation  time  r  =  fi/p  we  determine  the  binary-mixture  viscosity  using  the  Wilke  formula 
[3]: 


1  +  Gab 


Pb  Mg 

pa  Mb 


-1 


Pb  1  +  Gba 


pa  Mib 

Pb  M\.a 


-1 


where  Gab  = 


1  +  \Jpa/Pb\/Mb/Mg 

2^2  [1  + Mg/ Mb) 


(16) 


where  Mg  and  Mb  are  the  molar  masses  of  gases  a  and  b  respectively. 

Note,  that  the  collision  frequencies  and  the  binary-mixture  viscosity  are  parameters  “external”  to  the  QGDM 
model  and  could  be  determined  by  other  estimations. 


IV  SHOCK  WAVE  STRUCTURE  IN  A  HELIUM  -  XENON  MIXTURE 

As  a  first  example  of  using  the  QGDM  equations  we  considered  the  problem  of  a  stationary  shock  wave 
structure  in  a  mixture  of  helium  {He  is  gas  a)  and  xenon  (Xe  is  gas  b).  Density  profiles  for  these  gases, 
measured  with  the  use  of  an  electron  gun  and  a  laser  interferometer,  can  be  found  in  [11].  Measurements  were 
performed  for  the  following  variants: 

•  variant  VI  -  98,5%  He  and  1.5%  Xe 

•  variant  V2  -  97%  He  and  3%  Xe 

•  variant  V3  -  94%  He  and  6%  Xe 

•  variant  V4  -  91%  He  and  9%  Xe 

For  variant  V2,  there  is  a  calculation  by  the  DSMC  method  [1],  the  results  of  which  can  be  considered  as  a 
reference.  It  coincides  practically  with  the  experimental  data. 

The  system  of  equations  (7)  -  (9)  is  solved  in  non-dimensional  variables  taking  as  dimensional  scales  the 
following  characteristics  of  gas  a  in  the  upstream  flow:  pgre}  is  the  density,  acre/  =  -s/laRaTgref  is  the  sound 
velocity  at  temperature  Tgref,  Mref  is  the  mean  free  path,  that  is  computed  as  in  [1]: 
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FIGURE  1.  Density  (a)  and  temperature(b)  profiles  in  a  He-Xe  mixture 


ifj.  1 


(17) 


Then  relations  between  the  dimensional  and  dimensionless  parameters  have  the  following  forms  (all  param¬ 
eters  of  gas  6  are  scaled  by  parameters  of  gas  a): 

P  —  PParef  j  ®  —  ^^aref  i  ^  —  ^^aref  i  P  —  PParef  ^aref  ’  ^  ~  ^Pare  f  ^aref  ’ 


^are/ 

7a  77a 


TT, 


aref  i 


x\ 


aref  > 


t  =  t 


^are  f 
^are  f 


\3 

^are  f 


QGDM  equations  do  not  change  their  forms  after  the  process  of  scaling.  The  relations  between  the  parameters 
of  the  gases  (link  equations)  write  as 


aa  = 


fa 


laPa 

Pa 


rfi  laPb  Ra 

J-b  =  -^^-5- 

Ph 


Pa  =  T^\ 


-  Pbref 

Pb  —  - 

Paref 


t  aref 


Tbr 


e/ 


here  pbref  and  Tbref  are  the  viscosity  coefficient  and  the  corresponding  temperature  of  gas  h,  used  in  the 
viscosity-law  (15). 

Table  VI  presents  the  parameters  of  the  mixture  before  the  shock  wave,  chosen  in  accordance  with  the 
experimental  data  from  [11]  and  with  the  calculation  from  [1].  The  Prandtl  number  for  the  gases  is  constant 
and  equal  to  Pr  =  2/3. 

The  boundary  conditions  on  the  right  and  left  boundaries  were  taken  from  the  the  Rankine-Hugoniot  condi¬ 
tions  for  a  stationary  shock  wave  in  a  gas  mixture.  The  variables  on  the  right  of  the  discontinuity  are  computed 
as  follows 


P2  -  Pi 


U2  =  Ui 


(7-1-  l)Ma2 
2 -f  (7  -  l)Ma2  ’ 
2-f(7-  l)Ma? 


(7-1-  l)Ma2 


P2  =  Pi 


2'fMa?  —  7  -f  1 
7  -b  1 


(18) 
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FIGURE  2.  Mean  temperature  and  density  profiles  (a)  and  diffusion  velocities  (b)  in  a  He-Xe  mixture 


where  subscripts  1  and  2  refer  to  Rankine-Hugoniot  conditions  upstream  (1)  and  downstream  (2)  of  the  shock 
wave.  The  component  temperatures  are  found  from  the  equations  of  state. 

Assume  that  the  temperatures  and  velocities  of  the  components  before  and  after  the  shock  wave  are  equal,  and 
the  mass-fraction  of  the  components  after  the  transition  through  the  shock  are  unchanged  [na\/n\  = 
rn,\/n\  =  rn,2/n2)-  Thus  on  the  basis  of  conditions  (18),  the  parameters  of  each  component  of  the  mixture  are 
derived  from  the  ratios: 


Pa2  —  ^a^a2)  ^al  —  ^1  —  ^1 1  —  '^bl  — 

Pb2  =  mbnb2,  Ta2=Tb2  =  T2,  Ua2  =  Ub2  =  U2.  (19) 

The  initial  conditions  are  a  discontinuity  at  point  x  =  0: 

at  ^  0  Pa  —  Pal}  Pb  —  Pbl}  ^  —  T^li  —  Ub  —  Ui . 

at  a;  >  0  Pa  =  Pa2,  Pb  =  Pb2,  Ta  =  Tb=T2,  Ua  =  Ub  =  U2.  (20) 

To  solve  the  (7)  -  (9)  system,  an  explicit  difference  scheme  was  applied  where  the  steady-state  solution  was 
obtained  as  the  limit  of  a  time-evolving  process.  All  spatial  derivatives,  including  the  convective  terms,  were 
approximated  by  central  differences  (see,  for  instance,  [12]). 

The  problem  was  solved  using  a  uniform  spatial  grid  with  a  convergence  criterion  Cp^  =  10~®.  When  refining 
the  grid  by  a  factor  of  2,  the  differences  between  the  computational  results  were  extremely  small,  which  allows 
to  conclude  that  grid  convergence  has  been  reached. 

The  profiles  of  gas-dynamic  parameters  (velocity,  density,  temperature)  are  given  in  a  normalized  form  on 
the  basis  of  upstream  and  downstream  Rankine-Hugoniot  conditions.  In  this  case,  p  ^  {p  —  Pi)/{p2  ~  Pi)] 
similarly  for  the  temperature.  For  the  velocity  u  ^  (u  —  M2)/(mi  ~  M2)- 

Let  the  computational  results  for  variant  V2  be  considered  in  detail.  In  Figs.  1-3  the  profiles  of  gas-dynamic 
parameters  at  the  shock-wave  front  are  shown  in  comparison  with  the  corresponding  results  obtained  in  [1]  on 
the  basis  of  the  DSMC  method.  The  curves,  corresponding  to  the  DSMC  calculations,  are  superimposed  on 
the  QGDM  data,  so  that  the  values  of  the  mean  density  coincide  at  x  =  0. 

In  Fig.  1  are  the  profiles  of  density  and  temperature  of  helium  and  xenon.  In  Fig. 2a  are  the  distributions 
of  the  mean  density  and  temperature  for  the  mixture.  Like  in  the  DSMC  model,  the  mean  temperature  of  the 
mixture  is  very  close  to  that  of  helium;  and  the  temperature  of  xenon  overshoots  its  final  value  by  ~  10%  and 
the  mean  mixture  temperature  is  close  to  the  helium  temperature. 

In  Fig.  2b  are  presented  the  diffusion  velocities  Uda  and  Udb,  reduced  by  the  upstream  flow  velocity. 

Uda  =  Ua-U,  Udb  =  Ub  -  U,  (21) 
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FIGURE  3.  Diffusion  velocities  (a)  and  shock  wave  thickness  (b)  in  a  He-Xe  mixture 


where  u  is  the  mixture  velocity  that  is  defined  according  to  (13). 

In  Fig.  3a  the  xenon  concentration  is  presented.  Within  the  shock  wave,  it  falls  to  approximately  half  its 
initial  value. 

The  curves  demonstrate  that  the  QGDM  model  reflects  at  least  qualitatively  the  main  features  of  the  flow. 

Variants  VI, V3  and  V4  correspond  to  the  conditions  of  density  measurements  by  Walenta  [11].  The  compar¬ 
ison  of  experimental  results  with  those  obtained  by  the  present  numerical  work  leads  to  the  same  conclusion. 

On  the  basis  of  the  calculations  for  variants  VI  -  V4  the  shock- wave  thicknesses  Sne/^He  and  Sxel^Xe  have 
been  plotted  in  Fig. 3b  against  the  concentration  of  Xe  in  the  upstream  mixture  and  compared  with  the  results 
of  [11].  In  this  case,  the  shock- wave  thickness  is  calculated  as 

^  _  p2  -  p\ 
m&x[dp/ dx) 

The  mean  free  path  for  each  component  is  computed  according  to  (17)  for  parameters  of  each  gas  component 
before  the  shock  wave.  The  experimental  data  are  plotted  as  a  solid  line;  the  authors’  results  as  a  dashed  line. 
All  the  curves  are  represented  in  the  form  similar  to  [11].  For  variant  VI  (with  the  smallest  Xe  concentration), 
the  experimental  and  computational  results  coincide  practically.  With  increasing  Xe  concentration,  the  cal¬ 
culated  shock  wave  thickness  is  larger  than  the  experimental  values.  Nevertheless,  the  qualitative  behavior  is 
reproduced  by  the  calculations.  This  is  also  consistent  with  the  well-known  fact  that  the  relative  shock  wave 
thickness  increases  when  the  upstream  Mach  number  decreases. 

The  shock  wave  thickness  is  a  very  sensitive  characteristic  feature  of  the  problem,  and  its  calculation  based 
on  moment  equations  for  a  single-component  gas  corresponds  to  the  experimental  data  only  in  the  case  of  small 
Mach  numbers  Ma  <  2. 
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V  ARGON  -  HELIUM  DIFFUSION 

As  a  second  example  of  application  of  the  QGDM  equations,  the  problem  of  helium  and  argon  mass  diffusion 
was  studied  for  conditions  that  correspond  to  a  computation  by  the  DSMC  method  [1].  Let  two  reservoirs, 
filled  with  the  gases  be  located  at  a  distance  L  =  Im.  He  is  gas  a  in  the  right  reservoir,  and  Ar  is  gas  h  in 
the  left  reservoir.  The  number  densities  in  the  reservoirs  are  kept  constant  and  equal  to  n  =  2.8  •  10^°m“^.The 
gases  in  the  reservoirs  are  assumed  to  have  the  same  temperature  T  =  273K  and  the  same  velocity  equal  to 
zero. 

Using  these  constants,  the  missing  initial  data  can  be  obtained:  helium  density  pa  =  nnia  =  1.862  • 
10“®kg/m^;  sound  velocity  Oa  =  -\/7a RaTa  =  971.9m/s;  mean  free  path  computed  by  formula  (17)  Aa  = 
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FIGURE  4.  Number  densities  and  diffusion  velocities  in  the  Ar  -  He  diffusion  problem. 


1.479  •  10“^m;  argon  density  =  nrrib  =  1.856  •  10“®kg/m^;  sound  velocity  05  =  \/^jbRbTh  =  307.81m/s;  mean 
free  path  computed  by  formula  (17)  is  =  4.63  •  10“^m. 

As  in  the  previous  section,  computing  was  performed  in  dimensionless  variables  with  all  quantities  normalized 
by  the  parameters  of  gas  a  -  helium  in  the  reservoir. 

A  one-dimensional  plane  flow  described  by  equations  (7)  -  (9)  was  considered.  As  boundary  conditions  the 
following  non-dimensional  relations  were  used: 
at  the  left-hand  boundary  (1) 


^  _  1  _  iA— 10  _  _  1  _  p. 

Pa  —  J-*  J-U  ,  Pb  —  ,  la  —  lb  —  -!-•)  —  U, 

at  the  right-hand  boundary  (2) 

« =  1.-10-“  ,.  =  10-“.  r.  =  r.  =  i..  ^  =  ^  =  0. 

That  is,  we  assumed  that  in  each  reservoir  a  fraction  ~  10“'^°  of  the  other  gas  exists.  At  initial  time,  the 
density  of  the  components  between  the  reservoirs  is  assumed  to  change  linearly  : 


Pa{x)  =  ^  +  =  0), 

«(.)  =  >’'■(- =  ^> -»>(-  =  »>,  +  =  0). 

We  used  the  same  numerical  algorithm  as  in  the  previous  section  when  solving  the  QGDM  equations.  The 
problem  was  solved  using  a  uniform  space  grid  consisting  of  339  points  with  spatial  grid  step  h  =  0.2  ,  which 
corresponded  to  0.2Aa  and  0.64A;,. 

The  number  densities  and  diffusion  velocities  of  both  gases  are  plotted  in  Fig.  4  against  the  position  between 
the  reservoirs.  In  both  figures  the  comparison  with  DSMC  results  [1]  is  given.  Again  the  present  results  agree 
at  least  qualitatively  with  the  reference  results.  The  point  of  equal  concentration  point  is  shifted  from  the 
middle  of  the  domain  to  the  left,  closer  to  the  reservoir  containing  the  heavier  gas.  The  diffusion  velocity  of 
helium  is  larger  than  that  of  argon.  The  diffusion  velocity  of  helium  exhibits  a  minimum  in  the  middle  of  the 
computational  domain. 


303 


TABLE  1.  Dimensional  parameters  of  mixture  components  for  variants  VI  —  V4 


VI 

V2 

V3 

V4 

He  Xe 

He  Xe 

He  Xe 

He  Xe 

p(kg/mQ  •  10^ 

5.15  2.57 

5.16  2.22 

4.91  10.3 

4.57  14.8 

p(Pa) 

33.14  0.51 

33.21  1.02 

31.62  2.02 

29.42  2.91 

r(K) 

1  310. 

n(m/s) 

3076.76 

2882.6 

2672.8 

2530.3 

Ma 

2.97  17.01 

2.78  15.93 

2.58  14.78 

2.44  13.99 

VI  CONCLUSION 

The  macroscopic  system  of  the  QGDM  equations  is  constructed  on  the  basis  of  the  kinetic  equation  system 
in  the  relaxation  approximation  to  describe  the  non-reacting-gas  mixture  flow.  Contrary  to  some  widely  known 
models,  the  QGDM  system  consists  of  the  equations  for  density,  momentum  and  energy  of  each  component, 

i.e.  it  is  a  two-fluid  approximation,  that  makes  it  possible  to  describe  in  detail  the  behavior  of  each  component. 

The  momentum  and  energy  equations  include  exchange  terms,  that  allow  for  the  appropriate  exchanges 
between  the  gas  components.  To  calculate  these  terms  the  cross-collision  frequencies  must  be  estimated.  The 
viscosity  of  the  mixture  must  also  be  estimated.  The  QGDM  model  includes  diffusion  processes  but  does  not 
require  coefficients  of  thermo-,  baro-  and  self-  diffusion,  which  are  included  in  the  Navier-Stokes  models,  and 
the  determination  of  which  is  a  separate  task. 

The  numerical  algorithms  developed  on  the  basis  of  the  QGDM  model  appear  to  be  more  stable  than  similar 
algorithms  based  on  the  conventional  conservation  equations. 

Thus,  the  QGDM  equations  are  worth  being  studied  further  for  future  other  applications. 
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